ALTERNATIVE SOLUTIONS OF INHOMOGENEOUS SECOND ORDER LINEAR DYNAMIC 

EQUATIONS ON TIME SCALES 



DOUGLAS R. ANDERSON AND CHRISTOPHER C. TISDELL 

Abstract. We exhibit an alternative method for solving inhomogencous second-order linear ordinary dynamic equa- 
tions on time scales, based on reduction of order rather than variation of parameters. Our form extends recent 
(and long-standing) analysis on R to a new form for difference equations, quantum equations, and arbitrary dynamic 
equations on time scales. 



1. Introduction to second-order ordinary dynamic equations 

A very common equation in mathematics, mathematical physics, and engineering is the inhomogencous second- 
order hnear ordinary differential equation 

: (1.1) y"it)+pit)y'{t)+qit)yit)^rit), teR, 

its ordinary difference equation counterparts 

(1.2) A{Ay)it)+pit)Ayit)+qit)yit)=rit), t e Z, Ay{t) := y{t + I) - yit) 

or 

: (1.3) y{t + 2)+a{t)y{t + l) + l3{t)yit)=rit), t e Z, 

the related ordinary quantum equation 
■ (1.4) Dh{Dhym+p{t)Dhy{t) + q{t)y{t)=r{t), teh^, h>l, Dhy{t) 

or the recently introduced ordinary dynamic equation on time scales given by 

(1.5) y^'^(t)+p(%^(t)+q(%(t)=r(i), ^6T'^^ 

where the differential/shift operators in (jl.ip — (jl.Sp represent differentiation with respect to t on the corresponding 
time scales, respectively. In the general setting represented by (jl.Sp . the functions p, q, and r are real- valued 
right-dense continuous scalar functions of t satisfying the regressivity condition 

1 + /i(i) [~p{t) + il{t)q{t)] + 0. 
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Recall that on a time scale T, namely any nonempty closed subset of the real line, the delta derivative is given by 

s^t a{t) ~ s 

provided the limit exists, where a is the forward jump operator cr(t) :— mi{s G T : s > i}, and y"^ — y o a. The 
graininess fi is simply fj,{t) = (7{t) — t. For more on time scales see Hilger [7]. 

Indeed, extensive analysis of (jl.5l) and its solution forms can be found in Bohner and Peterson ^3^, while the 
ordinary differential equation (|l.ip is studied by, for example. Blest [2], Boyce and DiPrima [4], Hille [8], Ince [9], 
Johnson, Busawon, and Barbot [TT], Kelley and Peterson [T3], whereas the difference equation (|1.2p appears in 
Agarwal PP, Elaydi [5], and Kelley and Peterson [T^]. A commonly used technique to solve (jl.ip — (jl.5|) is Lagrange's 
variation of parameters method. In this approach, a solution y of (11.11) — (jl.Sp takes the form y = yu + yd, where yu is 
the complementary solution of the corresponding homogeneous or undriven (r = 0) form of (|l.ip — (|1.5p . and yd is any 
particular solution of the inhomogeneous or driven equations (jl.ip — (jl.Sp . If yi and y2 are two linearly independent 
solutions of the corresponding homogeneous equation, then it is well known that we may write y — cij/i + C2J/2 + yd 
for arbitrary constants ci and C2. For example, using variation of parameters, the form of a particular solution for 
(fTS]) is O Theorem 3.73] 



which reduces to 

J a W{yi,y2){s) J a W{yi,y2){s) 

for pTTj) . and to 

/,^ yi{s + l)r{s) y2{s + l)r{s) 

^M^(2/i,2/2)(s + l) ^ W^(2/i,y2)(s + 1) 

for (|1.2p . where in each case W{yi,y2) is the Wronskian of yi and y2, defined appropriately for each time scale. 
Notice that the integrals/summations always involve both yi and y2- Recently, Johnson, Busawon, and Barbot |11] 
derived two alternative solution forms for y^i for the ordinary differential equation (|l.ip . namely 



(1.7) Vdit) = y^{t) I -dt^ i e M, * = 1, 2. 



As the authors point out in [TT], either yi or 2/2 may be chosen in (11.71) . depending on which one yields an easier 
integral to compute. 

Unfortunately neither the technique nor the results in [11] are new, as they are discussed in Blest [2] and the 
undergraduate textbook by Boyce and DiPrima {4, Exercise 28, pl85], to cite just two examples. In fact, Yosida 
[l4l p. 29] calls the technique D'Alembert's reduction of order for ordinary differential equations, while according 
to Jahnke [TIT, p. 332], the method of reduction of order dates back at least as far as Euler (1750). In addition, 
Clairaut, Lagrange and Laplace have been involved with the method and thus it may be impossible to attribute to 
one mathematician [TU'j p. 332]. In any case, the method is at least 250 years old. 
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Despite the non- novelty of (I1.7P for the ode such a resuh would be new on general time scales for (11.51) . 

Consequently one of our goals in this paper is to generalize and extend (|1.7p to an alternative form for particular 
solutions to the ordinary dynamic equation on time scales ()1.5p . which would then nevertheless result in new alter- 
native forms for the difference equations (|1.2I) or (|1.3p . and the quantum equation (|1.4p as simple corollaries, as well 
as contain (|1.7p for (11.11) . This development is different from that given by Bohner and Peterson j3l Chapter 3]. 

An additional goal of this paper is to present some simple results that guarantee the nonmultiplicity of solutions 
to initial value problems associated with (jl.Sp . Our approach is based on simple inequalities and does not require a 
knowledge of matrix theory, in contrast to [3] Theorem 3.1; Corollary 5.90]. 

In Section [2] we address the question of nonmultiplicity of solutions, while in Section [3] we develop a method for 
solving inhomogeneous second-order linear ordinary dynamic equations on time scales, based on reduction of order 
rather than variation of parameters. Section contains some special cases that illustrate our results. 

2. NONMULTIPLICITY 

In this section we consider the notion of nonmultiplicity of solutions to the linear initial value problem (jl.5p with 
initial conditions 



where A,B(^M. and to G T. Let /CM and to E Ij := I T. By "nonmultiplicity of solutions" we mean that 
our theorems will present conditions under which (|1.5p . ()2.ip will have, at most, one solution y — y[t) for t > to, 
t^If. Such information is highly valuable, for example, when constructing explicit solutions to problems as we can 
determine when the constructed solution (or unique linear combination of solutions) will be the only solution to the 
problem at hand. Our techniques follow those of Coddington O Ch.2, Sec. 3]. 

In the following, we will make use of the time-scale exponential function defined in terms of right-dense continuous 
functions p satisfying the regressivity condition 1 + ix{t)p{t) ^ for all t G T''. Given such a p, the delta exponential 
function [3l Theorem 2.30] is given by 



where Log is the principal logarithm. It follows that ep(i, a) is the unique solution to the initial value problem 
0^(t) —p{t)<j){t), (l)[a) = 1 on T. We will denote l/ep(t, a) by eQp{t,a). 
We will require the following lemma from 3j Theorem 6.1, p. 255]. 

Lemma 2.1. Let £ e 71+ and v e Cjj(T), and let to G T. Ifv^{t) < i{t)v{t) for all t > to, then v{t) < v{to)ei(t,to) 
for all t > to- 



(2.1) 



y{to) = A, 



y'^ito)^B, 
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Our initial analysis will concern the case of (jl.Sp with constant coefficients, namely 

(2.2) y^^it)+py^it)+qyit)=rit), y{to) = A, y^{to) ^ B, 

where p,q gM. are constants. The following result involves the homogeneous form of (|2.2|) . that is 

(2.3) y^^{t)+py^{t) + qy{t)^0, 

and provides an estimate on the growth rate of solutions to (j2.3l) . For this estimate, we define 

\\y{t)h:={{y{t)r + {y'^{t))'y^\ t e h, t>to. 
Theorem 2.2. Let k -.^ 1 + \p\ + \q\. If y is any solution to (|2.3p on then 

(2.4) \\y{t)\\2 < \\y{to)hek{t, to), t e h, t> to- 
Proof. Let u{t) — \\y{t)\\2, where y is a solution to (|2.3p . For all i G we have 

u^{t) = (2/(t)+2/-(t))y^(t) + (y^(t)+y^-(t))y^^(t) 

= {2y{t) + nit)y^{t)) y^{t) + {2y^{t) + Kt)y^^it)) y^^(i) 
= 2yit)y^it) + ^i{t){y^{t)f + 2y^(i)y^^(<) + /.(t)(y^^(t))^ 

Now, apply Young's inequality 2ab < a? + to the first term and replace y^'^ with —py^ — qy to obtain 

«^(<) < mf + {y^ityf + M(t)(y^(t))^ + 2y'^{t)[~py^{t) - qy{t)] + M(i)hra^(t) - qvit)]' 

< (1 + 2|p| + \q\) {{y{t))' + (y^(t))2) + ^(<) [l +p2 + |p||^| + ^2] ((j^(^))2 ^ (yA(^))2) 

< 2(1 + bl + \q\) {{y{t)f + iy^{t)f) + /.(<) [l + 2(b| + |g|) +/ + 2|p||g| + q'] {{y{t)f + iy^{t)f) 
= {k®k)u{t). 

Above, © is known as the "circle plus" operator, z (B w := z + w + ^xzw, see [31 p. 54]. 
Thus, the conditions of Lemma I^TT] hold with v ~ u and £ — k k. Consequently, 

u(i) < u(to)efeefe(i,to) = M(to)(e/c(i,to))^, t>to, t e h, 

and therefore 

\\y{t)h<\\y{to)hek{t,to), t>to, teij. 

Theorem 12.21 now leads to the following nonmultiplicity result for solutions to (|2.2p , (|2.ip . 

Theorem 2.3. Let r G Q,-ca(Ij^ ;IR), cmd t^ G /j. T/ie dynamic initial value problem (|2.2p . (|2.ip /las, at most, 
solution y — y(t) for t > to with t £ Ij. 
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Proof. Assume there are two solutions x and y. Let z{t) :~ x(t) — y{t) for t e Ij, and note that z must satisfy 
the homogeneous equation p.3p together with the homogeneous initial conditions z{to) = 0, z^{to) — 0. Now, by 
Theorem 12.21 we have 

\\z{t)\\2<\\z{to)\\2ekit,to)^0, t>to, teh, 
and thus z{t) =0 for alH > to, ^ £ -^t- It follows that x — y. □ 

The following result is an extension of Theorem 12. 21 and concerns estimates on solutions to the homogeneous form 
of (|1.5p with variable coefficients, namely 

(2.5) y^^{t)+p{t)y^{t) + q{t)y{t)=0, t e 

Theorem 2.4. Let t £ Ij. Let pi,qi e M 6e nonnegative constants such that 

\p{t)\<Pl, \q{t)\<qi, telf, t>to, 
and let ki := 1 +pi + qi. If y is any solution to (|2.5p on If, then 

Mt)\\2<\\y{to)hekAiM). t>to, t€h. 

Proof. Our proof follows similar lines as that of the proof of Theorem l2.21 and thus we just sketch the details. Letting 
u{t) — \\y{t)\\2, we obtain 

u'^it) < {y{t)f + {y^{t)f + fi{t){y^{t)f + 2y^(t)[-p(t)y^(t) - q{t)y{t)] + W - q{t)y{t)f 

< 2(1 + pi + qi)u{t) + ^l{t) [l + 2{pi + qi) + pI + 2piqi + qj] u{t) 

< {ki®ki)u{t). 
Thus, applying Lemma |2. II we obtain 

w(i) < ^t(^o)(efcl(^,^o))^ * > ^o, t € h, 

and therefore 

\\y{t)\\2<Mto)\\2ekAtM). t>to, teh. 

□ 

Theorem 12.41 now leads to the following nonmultiplicity result for solutions to (|1.5|) . (|2.ip . 

Theorem 2.5. Lett G Ij. The dynamic equation (|1.5p with initial conditions (j2.ip has at most one solution y — y{t) 
for t > to with t Cz Ij. 

Proof. Assume there are two solutions x and y, and let z{t) = x{t) — y{t) for t > to with t G It- Note that z 
satisfies (|2.5p and the initial homogeneous initial conditions z{to) ~ 0, z^{to) = 0. Since Ij may be unbounded, the 
coefficient functions p and q may not be bounded on Ij, and so Theorem 12.41 mav not be directly applied to z. We 
let t be any point in Ij such that t > to, and let Jj be any closed, bounded interval of Ij such that Jf has to as a 



6 ANDERSON AND TISDELL 

left cndpoint and Jt contains t. On Jt, p and q are both bounded, say by pi and qi, respectively. We can now apply 
Theorem 12.41 to z on Jt, and so z{t) — for all t E Jj. Now since t was chosen to be any point in Ij with t > tg, we 
have shown that x{t) — y{t) for all t E Ij with t > tg. □ 

Remark 2.6. The quest for nonmultiplicity of solutions on intervals to the left of to is a more delicate affair. The 
results of this section may be extended to include nonmultiplicity of solution for t < by adapting the proofs and 
obtaining inequalities like 

u^{t) > ~2ku(t), for all t<tQ, t€ h. 

However, there is a price to pay - the graininess function of the time scale would need to be bounded above. This is 
due to regressivity coming into play. 

3. ALTERNATIVE SOLUTION FORMS FOR THE GENERAL INHOMOGENEOUS DYNAMIC EQUATION 



In this section we state and prove the main result, namely a new form for a solution of (jl.Sp in the spirit of (|1.7p . 
As mentioned previously, this will then give us a new form for the difference equation p.2p and quantum equation 
(fOl) as well. 



Theorem 3.1. Let T be a time scale and let p, q and r be real-valued right-dense continuous scalar functions of 
t G T with p and q satisfying the regressivity condition 

(3.1) i + fi{t)[-p{t) + f4t)q{t)]^0, teT\ 

For all t eT, let yi and y2 satisfy 

(3.2) yf^(t)+p(t)yf(i)+g(t)y,(<)=0, * = 1,2, 

and W{t) := yi{t)y2 {t) — y2{t)yf' {t) ^ 0. The general solution of the linear inhomogeneous second-order ordinary 
dynamic equation 

y^^{t)+p{t)y^{t)+q{t)y{t) = r(t), t 

is 

(3-3) y{t) = ciyi{t) + 02^2 ^) + yi{t) / ^ TTwTTn 

or 

(3-4 y{t) ^ ciyi{t) + C2y2{t) + y2{t) / ^ TTwTT^ 

where ci and C2 are arbitrary constants and e(_p_|_^g') (•, a) is the time-scale exponential function. 

Proof. Clearly p.3p and p.4p are of the expected form y = cij/i + C2y2 +yd- One could easily verify that a function 
of the form 

(3.5) y,{t) = y^it) [ it)eei-,^,,) i-jt) , «) At 



INHOMOGENEOUS SECOND-ORDER LINEAR EQUATIONS 7 

is a particular solution of (jl.Sp using the time scale calculus, but that would not give any insight into where (j3.5p 
comes from. Thus, to derive p.5|) . assume we have a particular solution to the inhomogeneous equation (|1.5|) of the 
form Udit) = yi{t)v{t), where yi solves the homogeneous equation p.2p and u is a function to be determined. Then, 
using the product rule {fg)'^ = gf^ + f^g^, we have 

yf{t) = v{t)yHt) + ynt)v^{t) 

= v{t)yHt)+y,{t)y'[{t)v^{t)/y,{t), 

and using the product rule again with the quotient rule (^^^ — ^~/^ we see that 

y, it) = v{t)y, {t)+y^ {t)v {t) + -^-^^ 

= v{t)y^^it) + yt%t)v^{t) + {y,it)y'[it)v^{t)f /yUt) - v^it)y^it). 
Since we are assuming yd is a particular solution of the inhomogeneous equation (jl.Sp , we must have 
r{t) = y^^{t)+p{t)yf{t)+q{t)y4t) 

= v{t) (yf^(t) +p{t)yHt)+q{t)y,{t)) + p{t)y,{M {t)v^ {t) / y,{t) 
+yt%t)v'^{t)+[y^{M{t)v^{t)\'' /yUt)-v'^{t)yHt). 
Now yi is a solution of the homogeneous equation (13.21) . so after simplifying we multiply by yf to get 

riMit) = yKt)p{t) [y^{M{t)v'^{t)\ /y^{t) + [y,{M{t)v'^{t)\'' + ynt)v^it) - (0) • 

Make the substitution u ~ yiyfv'^, then use the simple formula — f = fJ-f^ to get 

rmH) = ynt)p{t)u{t)Mt) + u^{t) + yf{t)v^{tMt)y^^{t) 

yi{t) 

Using the simple formula — f = nf^ again and rearranging, we see that 

u^{t) + {p{t)-f,{t)q{t))u{t) = r{t)yUt). 
Focusing on the coefficient of m, we note that 

p{t) - ^lit)qit) = -{-p{t) + p{t)q(t)) 

= - e {Qi^p{t) + nit)q{t))) 

e{-pit) + Ai(t)q(t))ee(_p+^q) {t, a) 



et(_,+^,)(t,a) 
^e(-p+p<?)(^'") 

Consequently we have that 



[1 + M(t) (e(-p(t) + ^i{t)q{tm ee(_p+^,) [t, a) 
°e(-p+/i<?)( 



(ee(--p+/<9)(i>a)w(i))'^ = r{t)yf{t)e1^(_p^^^){t,a). 
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Solving while recalling that u — yiy"v^ ^ we arrive at 

so that via yd — yiV we obtain p.Sp . □ 

Remark 3.2. The regressivity assumption in p.ip is not at all unusual, as it is automatic in the case T = M since 
/i = 0, and it is assumed in the variation of parameters theorem on general time scales; see [3j Definition 3.3]. The 
following theorem is a simple corollary of Theorem 13.11 

Theorem 3.3 (Reduction of Order). Assume p and q are real-valued right-dense continuous scalar functions of 
t G T satisfying the regressivity condition 

(3.6) 1 + fx{t) [-p{t) + fi{t)qit)] ^0, te T«. 

If yi is a solution of the linear homogeneous second-order ordinary dynamic equation p.5|) . then 

is a second linearly independent solution of (|2.5|) . Similarly, if y2 is a solution of (|2.5p . then 



(3.8) yiit)=y.{t) I fizig^2l^Ai 



is a second linearly independent solution of (|2.5p . 

Proof. We will prove p.7p . since the proof of (13.81) is similar. Thus, assume yi is a solution of (12. Sp . Since r{t) = 
in this case, and general antiderivatives are used in Theorem 13.11 we may choose the constant of integration in (|3.5p 
in such a way that 

J r{t)yiit)ee(^^p+f,g){a{t),a)At = 1, 

so that (13.51) becomes 



in other words a particular solution of (|2.5p . One could also verify p.7p directly using the delta derivative rules. To 
show linear independence, we calculate the Wronskian of j/i and y2, namely 

(3.9) W{yi,y2m = yi{t)y^ [t) - yt {t)y2{t) = e(_p+^,) (t, a) ^ 

for alH G T by the regressivity assumption p.6p : see |31 Theorem 2.48]. □ 

Remark 3.4. It is also possible to show the equivalence of the result shown in p.Sp with that due to the variation 
of parameters formula in equation (|1.6p . Let yi be a solution of the linear homogeneous equation (12. 5p . Assuming 
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p.6|) . and appropriating the Wronskian of yi and the form of 2/2 given in p.7|) as calculated in (\3.9\i . the solution 
given in ()1.6p is 

. ^ /■* yUs)r(s) , , , /■* 2;2('S)r(s) 



/q e(-p+^g)(cr(s),a) e(_p+^g)(cr(s),a) 



(3.10) -yi{t) / r(s)y^(s)ee(-p+^,)(a(.),a) / ^'^t^-/ A^A.. 

If we use the integration by parts formula / fg^At — f{t)g{t) — / f^g"'At, see O Theorem 1.77(vi)], on the integral 
in l|3.5l) . where we have taken 

fit)= fr{s)yUs)eei-p+,,M^),a)As and gA(^) ^ ^(^p±mM 



we get IPH})) . 



The next corollary concerns another possible second-order linear dynamic equation discussed by Bohner and 
Peterson [3, (3.6)]. 

Corollary 3.5. Let a, (3 and r he real-valued right-dense continuous scalar Junctions on T with a satisfying the 
regressivity condition 1 + fi(t)a(t) ^ for t eT"^. For all t eT, let yi and y2 satisfy 

(3.11) y^^it) + ait)y^^t) + mvUt) = 0, » = 1,2, 

with W{t) :— yi{t)y2 {t) — y2{t)y^{t) ^ 0. The general solution of the linear inhomogeneous second-order ordinary 
dynamic equation 

(3.12) y^^t) + a{t)y^%t) + mv^it) = r{t), t e T'^\ 



.„ioN f,^ f,^^ u^ I ee^{t,a) J r{t)y1{t)e^{a{t),a)At 
(3-13) y[t) ^ ciyi[t) + C2y2[t) + yi(t) / t—— At 



CiiA^ t,\^ , ee„{t,a) J r{t)yUt)ea{a{t),a)At 
(3.14) y{t) = ciyiit) + C2y2{t) + y2[t) / At, 



y2{t)y^(t) 

where ci and C2 are arbitrary constants, and where ea(-,a) is the time-scale exponential function. 

Proof. Rewriting p.l2p using the simple formula = f + /i/^, we see that we arrive at an equation of the form 
(jl.Sp . where 

p(t)^4^1±#l^ and m- ^^'^ 



l + ^i{t)a{t) ' ' 1 + fj.{t)a{t)' 

Then we see that the term — p + /ig in Theorem 13.11 above is given by 

~p{t) + ii{t)q[t) = - Qa{t), 

1 + fi{tja(t) 
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and this corollary follows. □ 

4. SECOND-ORDER LINEAR ORDINARY DIFFERENCE EQUATIONS 

In this section we recount the results of the previous section when T = Z, that is to say for second-order linear 
ordinary difference equations. The following corollaries to Theorem 13.11 are obtained by simply taking T = Z in 
Theorem 13.11 above. The first uses the forward difference operator form of the second-order linear equation ()1.2[) . 
and the second uses the shift form (11.31). 



Theorem 4.1. Letp, q andr be real-valued scalar functions oftEZ, withp and q satisfying the regressivity condition 
(4.1) 1 - p{t) + q{t) ^ 0, teZ. 

Let yi and y2 satisfy 

A2y,(t) + p{t)Ay,{t) + q{t)y,{t) =0, i = 1, 2, 

and W{t) := yi{t)Ay2{t) — y2{t)Ayi{t) ^ 0. The general solution of the linear inhomogeneous second-order ordinary 
difference equation 

A^y{t) +p{t)Ay{t) + q{t)y{t) = r{t), t G Z, 

is given by 



(4.2) y{t) = ciyi{t) + C2y2{t) + yi(t) ^ ■ 



(4.3) v{t)^c^yi{t) + c2y2{t) + y2{t)Y. 

and where ci and C2 are arbitrary constants. 



yi{t)yi{t + l) 



n;iI(i-p(j) + g(.))E nCa-^l^a)) 



Theorem 4.2. Let a, /? and r be real-valued scalar functions of t 1, with regressivity condition /3{t) ^ for t G Z. 
Let yi and y2 satisfy 

y,{t + 2) + a{t)y,{t + l) + p{t)y,{t)^0, i = 1,2, 

with W{t) := yi{t)y2{t + 1) — y2{t)yi{t + 1) 7^ 0. The general solution of the linear inhomogeneous second-order 
ordinary difference equation 

(4.4) yit + 2) + ait)yit + l)+(3it)yit)^r{t), i G Z, 

is given by 



(4.5) yit) = ci2/i(<) + C2j/2(0 + yiit) J2 ■ 



yi{t)yi{t + i) 



(4.6) yit) = ciyiit) + C2y2{t) + y2{t)J2 u^u 

^ y2[t)y2[t + l) 
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where Ci and C2 are arbitrary constants. 

In [TTl Section 3] the authors apphed their alternative formula (|1.7p to give a general treatment of (|l.ip in the 
case of constant coefficient functions p{t) = 2P and q{t) = Q with P and Q real constants with Q / and ^ Q. 
We follow suit by giving a general treatment of ()4.4|) with constant coefficients using the results above from Theorem 
14.21 To wit, we analyze the ordinary difference equation 

(4.7) y{t + 2) + 2ay{t + l)+ py{t) ^r{t), t G Z, 

where a and /3 are real constants with /3 7^ and ^ fi. The solution of (I4.7|) can be written as 

From Theorem 14.21 we know that a particular solution to the inhomogeneous equation (j4.7p has the form 



which simplifies to 



a:— a s—a 

Thus, complete expressions for the solution of (|4.7p are given by 

t-i x-i 

y{t) = ciA* + C2A* r{s)\\+^-^'^ P--'-^ 



X—a S—a 



y(t) = c,\{ + + ^ ^ r(s)A*+^-2-r-i-^ 



x—a s—a 



where for simplicity we have taken Ai — —a — ^ — fi and A2 ~ —a + yc?— /3. 

Remark 4.3. A general treatment of (jl.Sp with constant coefficients on arbitrary time scales and even quantum 
equations is made difficult by the fact that, as seen in (j3.5p . even with p and q taken to be constant functions, the 
term — p + iiq is not constant except for the very special cases of T = R and T = /iZ. 
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